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ABSTRACT 

Massive galaxies today typically are not forming stars despite being surrounded by 
hot gaseous halos with short central cooling times. This likely owes to some form of 
"quenching feedback" such as merger-driven quasar activity or radio jets emerging 
from central black holes. Here we implement heuristic prescriptions for these phe- 
nomena on-the-fly within cosmological hydrodynamic simulations. We constrain them 
by comparing to observed luminosity functions and color-magnitude diagrams from 
SDSS. We find that quenching from mergers alone does not produce a realistic red 
sequence, because 1 — 2 Gyr after a merger the remnant accretes new fuel and star for- 
mation reignites. In contrast, quenching by continuously adding thermal energy to hot 
gaseous halos quantitatively matches the red galaxy luminosity function and produces 
a reasonable red sequence. Small discrepancies remain - a shallow red sequence slope 
suggests that our models underestimate metal production or retention in massive red 
galaxies, while a deficit of massive blue galaxies may reflect the fact that observed 
heating is intermittent rather than continuous. Overall, injection of energy into hot 
halo gas appears to be a necessary and sufficient condition to broadly produce red and 
dead massive galaxies as observed. 
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1 INTRODUCTION 

Observations show a distinct bimodality in galaxy colours 
JStrateva et al.l|200ll;lBaldrv et al1l2004l;lBalogh et al.ll2004l ; 
iBell et all 12004 ; IWeiner et all 120051 ; IWillmer et all 12006ft . 

with massive galaxies today generally being ellipticals hav- 
ing little star formation or cold gas, while less massive galax- 
ies tend to be blue spirals with cold gas. Despite being rec- 
ognized in the earliest observations of galaxies, the origin 
of this bimodality remains poorly understood. The lack of 
galaxies in the region between the red sequence of ellipticals 
and blue cloud of spirals (often called the green valley) indi- 
cates that star formation must be quenche d rapidly in order 
to transform blue galaxies into red ones (|Bell et al.l [20041 : 
Blanton 2006). The preponderance of active galactic nuclei 
(AGN) in the green valley suggests that AGN are some- 
how connected to the proce ss of quenching star formation 
(e.g. [Schawins ki et all 120091 ). Several physical mechanisms 
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have been proposed to effect this transformation, with two 
frontrunners garnering much attention recently: feedback as- 
sociated with galaxy mergers, and quenching due to virial 
shock heating of accreted gas. Both of these invoke AGN as 
a crucial energy source, namely quasars in the former case 
and radio jets in the latter. 

Understanding the origin of massive, passive "red 
and dead" galaxies remains a problem for models of 
galaxy evolution. Hydrodynamic simulations of gas-rich 
galaxy mergers have been shown to induce starbursts and 
fuel p o werful quasars |Hernouisj [l989: Bar nes fc Hernouistl 
119921 : |Pi Matteo et al.l l2005l v Observations lend sup- 
port to the ide a of a connection be t ween mergers and 
powe rf ul AGN jHutchings et all ll98Sl; ISanders fe MirabeJ 



1996; 



Canalizo fc Stockton! 1200 ll : iBennert et all |2008| : 
but not all: e.g. iBahcall et al.l Il99 



Urrutia et al.l 200^. . _ 

McLure et ai]ll999l: iGrogin et all 12003 : iGabor et all 1200 



Georgakakis et al.l 120091 ') . In merger simulations, the com- 



bination of feedback from star formation and black hole 
accretion unbinds the gas from the merger remnant, leav- 
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ing a red elliptical ga l axy devoid of fuel for star forma- 
tion JSilk fc Reed 1 19981 ; ISpringel et aill2005l ; iHopkins et ail 
120081 ). This is now colloquially known as "quasar mode" 
feedback. Observations support various aspects of this pic- 
ture, incl uding high-velocity o utflows from post-starburst 
galax ies (|Tremonti et al.1 120071 ) and QSOs ijFeruglio et al.1 
120101 ). and a peak in the fraction of AGN host galaxies 



Silverman et al. 


2008 
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Schawinski et al 1 I2009I . but see also 
Xue et al.ll2010l ). Recent work hints 



that secular processes may induce AGN (|Cisternas et al.l 
but it is not clear whether such AGN alone can 
quench star-formation in their host galaxies. 

The presence of hot X-ray emitting gaseous halos 
around red and dead galaxies in groups and clusters sug- 
gests that these may also be associated with quenching. 
In this scenario, gas falling into the halo from the inter- 
galactic medium (IGM) collides with stationary hot gas al- 
ready in the halo. The resulting shock heats the gas to the 
virial temperature, convertin g the gravitationa l energy of 
infall to thermal energy (cf. feggen et al.l|l962l ; IsIIkl Il977l ; 
iRees fc Ostrikerlll977l ; Twhite fc Reeslll97sl y Several recent 
hydrodynamic simulations show that radiative cooling dom- 
inates over this shock heating in dark matter haloes below 
a rough critical mass of ~ 1O 12 M , which is suggestively 
close to the mass at which the color bimodality divides 
galaxies. Above this critical mass, a stab le shock forms and 
the halo develops a stable ho t gas halo l|Birnboim fc Dekel 
l2003l ; lKeres et al.ll2005l , l2009bl '). This hot halo will shock heat 
newly accr eting gas, thus stalling the ultimate fuel for sta r 
formation |Dekel fc Birnboiml [200r3 ; iBirnboim et al.ll2007h . 
Dense filaments of gas may penetrate this hot halo, avoid- 
ing the virial heating and providing fuel for star formation, 
bu t this is predicted to be uncommon at redsh i fts less than 
2 (iDekel fc Birnboiml l200rj ; lOcvirk et al. I l2008l ; IDekel" et al.l 
l2009al ) because filaments are less dense at late times. 

In massive galaxy cluster haloes w ith peaked densit y 
distributions (roughly half of all clusters; iBauer et al.ll2005t ). 
the hot gas appears to be cooling rapidly via X-ray emission 
I Fabian et al.|[l984r ). yet UV, optical, and NIR observations 
indic at e only trace levels of s tar formation ijSmith et al.l 
19971; [Crawford et al l 1 19991; iHicks fc Mushotzkvl 120051 ; 
Guillen et al.l 120081; iDonahue et alj 120101 1 and little cold 



gas (e.g. ISalome et al.l 120061 ') . Hence some additional heat- 
ing mechanism must counteract the X -ray cooling; this is 
the well-known cooling flow problem (|Fabianl [l994'). Ob- 
servations of X-ray cavities thought to result from pow- 
erful radio jets have motivated a picture where a ra- 
dio AGN in the central cluster gala x y provides this ex- 
tra heating dMcNamara et al.1 |2000| ; iFabian et al] l200d ; 
iMcNamara et all 120051 ). As gas cools from the hot halo, 
some of it accre tes onto the massive cen tral black hole, in- 
ducing an AGN (|Ciotti fc Ostrikerll 19971 ). The kinetic power 
emer ging via radio jets creates bubbles of typical size ~ 10 
kpc (jBirzan et al.l 120041 ') in the central regions of the clus- 
ter. These bubbles then expand and degenerate into sound 
waves that effectively isotropize outflow ene rgy and heat 
the hot intracluster gas on scales >30 kpc (IFabian et al. 
20031 ; IVoit fc Donahue! 120051 ; IFabian et alTl2006l ; iDe Young 
20101 ). From gas cooling to gas heating by the AGN, this 
entire cycle may repeat itself on ~ 10 8-9 yr timescales 



jCiotti fc Ostrikerll200lf). This is now referred to as "radio 
mode" feedback (|Croton et al.ll2006l ). 

Mechanisms unrelated to supermassive black holes may 
provide the additional heating in hot X-ray haloes. Ther- 
mal conduction in clusters provides central heating and in- 
hibits cooling instabilities, altho ugh it is not sufficient to pre- 
vent cooling flows in all clusters jZa kamska fc Naravanl2003l ; 
IParrish et al.ll2009l ; lRuszkowski et alj|2010h . Infalling galax- 
ies or stellar clumps may dynamically he at haloes by trans- 
ferring gravitational energy to the gas (IDekel fc Birnboiml 



20081; iKhochfar fc Ostrikerl |200S|; Ijohansson et al.l |2009[ : 



Ruszkowski fc Ohl2010] ; lBirnboim fc Dekel2010f l. Whatever 
the mechanism, the key requirement is to heat up the cores 
enough to prevent substantial cooling and star formation, 
while allowing central cooling times less than the Hubble 
time (or even < 1 Gyr) to persist in a substantial fraction 
of clusters and groups. 

Debate persists about the effectiveness of these pro- 
posed quenching mechanisms on scales smaller than a 
kpc, within a single galaxy. Various simulations use 
different models for black hole accretion and feedback 
dSpringel et~aH 120051; iBooth fc Schavd 120091 ; iDeBuhr et ail 
120091 ; iLevine et al.1 |201oT T and recent work raises doubts 
that even a powerful quasar cou ld significantly affect star - 
forming gas throughout a galaxy (|DeBuhr et al.ll2"009l . l2010T ). 
In radio mode feedback, the observed AGN energy out- 
put is generally sufficien t to counteract the radiative cool- 
ing in clusters or groups l|McNamara et al. I l2006l ; iBest et al.l 
120061 ; iGiodini et alJlioibl ). How this energy is efficiently dis- 
tributed over space and time to prevent substantial cooling 
and star formation remains unclear, and cosmological mod- 
els generally do not include realistic AGN radio jet physics. 
Hence while supermassive black holes are likely to be re- 
sponsible for quenching, direct ab initio simulations of this 
process on a cosmological scale remain beyond current reach. 

Semi-analytic models (SAMs), which marry analytic 
prescriptions for baryon physics to merger trees from cos- 
mological simulations of dark matter, have successfully in- 
corporated quenching mechanisms to create reasonable red 
sequence populations. Among analytic prescriptions for gas 
cooling and star formation, these models typically incorpo- 
rate a starburst and bulge growth during galaxy mergers. 
These merger prescriptions are not enough to make mas- 
sive galaxies red, so modelers include some suppression of 
cooling in hot , mass i ve haloes, usually v i a AGN feedback 
(ICroton et al. I 120061 : iBower et al. I l200rj ; ISomerville et al.l 
120081 ). While most SAMs contain such prescriptions, the ex- 
act amount of quenching feedback required and the mech- 
anisms by which they operate differ widely among models, 
possibly because these models must make many assump- 
tions regarding the gas dynamics which c an be quite uncer- 
tain (|Lu et al.ll2O10l ; lDe Lucia et al.ll2010l ') . 

Hydrodynamic cosmological models follow the gas dy- 
namics directly, but without some explicit quenching mech- 
anisms they generally do not yield a realistic red sequence 
with appropriate luminosity or mass functions. These mod- 
els greatly reduce the number of free parameters relative 
to SAMs by tracking the dynamics of gas inflows and out- 
flows directly, at large computational cost. Although hy- 
drodynamic models directly track the dynamics of galaxy 
mergers and yield hot gas in massive haloes, they still pro- 
duce only constantly growing, blue, star-forming galaxies at 
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the massive end of the mass function. Recent work has be- 
gun to include black hole fueling and feedback. These works 
have focused on reproducing the black hole mass-bulge ve- 
locity dispersion (M — a) relation, black hole properties, 
and AGN population properties, with significant successes 
fe.g ISiiacki et al]l2007f : IPegraf et ai1l201ol : iBooth fc Schavd 
120091 . 120101 )7 Typical models involve Bondi-Hoyle accretion 
of surrounding gas by the central black hole, with some as- 
sumed fraction of the accreted rest-mass converted to heat in 
the s urrounding gas (|Springel et al.l 120051 ; IBooth fc Schavd 
2009). These prescriptions, while pioneering, are poorly con- 
strained and may n ot properly model the dominant ac- 
cretio n mechanisms |Booth fc Schavel 120091 ; IPeBuhr et al.l 
2009). They also do not focus on producing red sequence 
galaxies as observed, and have had only limited success do- 
ing so. 

In this work we also employ hydrodynamic simulations, 
but we focus on reproducing the observed galaxy population 
rather than the observed black hole and AGN population. 
Our approach is therefore somewhat different, in a sense in- 
termediate between fully hydrodynamical simulations and 
SAMs. Here we simply ask, how is quenching related to 
the evolutionary properties of galaxies and their halos? We 
apply heuristic models for quenching by adding energy in 
various forms to galaxies and their surrounding gas during 
the evolution of the simulation, and ask which models are 
successful at producing red and dead galaxies as observed. 
We do not try to ex plicitly grow b lack h oles and account 
for their feedback. In IGabor et~ail (|2010T ). we conducted a 
precursor study on post-processed star formation histories, 
and concluded that both radio mode and quasar mode feed- 
back could in principle produce red galaxies if all future star 
formation was assumed to be suppressed, but several other 
attempted mechanisms such as shutting off hot mode accre- 
tion or recycled wind mode accretion could not. As such, we 
now focus particularly on the two popular mechanisms of 
radio mode and quasar mode, with the goal of determining 
which one drives the quenching of massive galaxies. 

Here, as in IGabor et all (|201Ch . we compare our sim- 
ulation results to observations of the galaxy population at 
redshift z ~ 0, and leave detailed evolutionary studies for 
the future. We consider the red sequence in color-magnitude 
diagrams, and the corresponding red galaxy luminosity func- 
tion. Due to uncertainties in dust models, we use galaxy 
stellar mass functions to compare our simulated blue cloud 
population observations. Our main result is that superwinds 
induced by galaxy mergers cannot create a red sequence, be- 
cause galaxies continue to accrete gas from the IGM even 
after a merger, while our quenching mechanism based on 
hot halo gas produces a successful match to the observed red 
galaxy luminosity function. This favors radio mode feedback 
as the primary driver for the formation of the massive red 
galaxy population. 

We describe the physics of our simulations in iJU includ- 
ing our newly implemented quenching models in tl2.il In SJ3] 
we describe the observational data. We compare the results 
of these new models to observed color-magnitude diagrams 
(CMDs) and luminosity functions (LFs) for local galaxies in 
<J4] We then consider physical consequences of our models 
in 5j5] and discuss future directions in the discussion. 



2 SIMULATIONS 

Our simulations are run with an extended version of the 
N-body + smoothed particle hydrody namics (SPH) cos- 
mological si mulation code GADGET-2 (|Springelll2005l ). as 
described in lOppenheimer fc Pavel (|2008T ) . Along with gas 
cooling and star-formation, our version includes galactic 
winds tied to the star formation and a chemodynamical 
model with sources including AGB stars and Type la super- 
novae. In the following subsections, we describe additional 
modifications to the code to implement quenching of star 
formation. 

For star formation in cold, dense regions, w e use 
the two-phase m odel of [S pringcl & Hcrnquist] 1)20031 ). This 
model, based on iMcKee fc Ostrikerj (|197tI v" assumes that 
any gas particle dense enough to be Jeans unstable con- 
tains a hot medium and cold clouds in pressure equilibrium. 
Stars form from the the cold clouds with a characteristic 
timescale, r. This ti mescale effectively sets the e fficiency of 
star formation, and ISpringel fc Hernquistl (2003) find that 
t ~ 2 Gyrs yields galaxy models t hat match the observed 
Kennicutt relation (jKennicuttll 19981 ) . The three model com- 
ponents exchange mass, energy, and metals via condensation 
and feedback from supernovae. 

Star-forming gas particles in the simulations spawn star 
particles with a probability based on the star-formation rate. 
The resulting collisionless star particle takes on a mass equal 
to half the original mass of a gas particle. The average gas 
particle spawns up to two star particles in this stochastic 
way. 

Along with the feedback implicit in the star-formation 
model, our simulati ons include an explicit stoch astic model 
of galactic winds (|Oppenheimer fc Pavel |200S| ). We eject 
star-forming gas particles from galaxies by giving them a 
velocity kick v w (|Springel fc H crnquist 2003). The probabil- 
ity that a gas particle will enter a wind in the first place is 
the mass-loading factor r\ (the ratio of mass entering winds 
to the mass of stars formed) times the probability for star 
formation. In our simulations v w is proportional to, and r\ 
inversely proportional to, the galaxy circular velocity. These 
scalings, whic h arise in the momentum-driven wind model of 
IMurrav et al.l |2005|), are motivated by observation s of winds 
in local galaxies " (|Martinll2005l ; IRupke et al.l 120051 ) . Once a 
gas particle enters a wind, we decouple it hydrodynamically 
until it reaches a density 0.1 times the critical density for 
star formation (up to a maximum duration of 20 kpc/wu,). 
This decoupling is intended to represent chimneys through 
which material can escape easily but that cannot be reliably 
modeled by the spherically-averaging SPH algorithm. 

In order to study the chemical enrichment of galaxies 
and the IGM, our code includes mod els for stellar mass loss 
and T ype la supernovae, described in lOppenheimer fc Pavel 
(2008). As the stellar population represented by a star 
particle ages, mass loss from AGB stars injects metal- 
enriched mass to the nearest 3 gas parti c les ba sed on the 
mass-loss rates from iBruzual fc Charlotl (|2003t). We esti- 
mate Type la supernov a rates for each star particle using 
IScannapieco fc Bildstenl (|2005l ). with a prompt and delayed 
component. Each supernova injects both energy and metal 
mass into the surrounding gas. 

With this treatment of galactic winds, our simula- 
tions match broad observations of star-forming galaxies and 
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the IGM, i ncluding the chemical enrichment of the IGM 
at z > 2 (|Oppenheimer fc Pavel 120061. 120081), the galaxy 
mass -metallicitv relation l|Pave et al .1120071 : iFinlator fc Pavel 
2008) , and the luminosity functions of h igh-redshift galaxies 
I Pave et all 120061 : IFinlator et al.1 120071 '). These simulations 
do not, however, yield massive red and dead galaxies, and 
galaxies are far too massive and bright at low redshifts. We 
require some additional mechanisms to halt star-formation 
in massive galaxies - the hydrodynamics and feedback mod- 
els we have described above are insufficient. 



2.1 Quenching Prescriptions 

We incorporate two independent quenching models (which 
may be combined) into GAPGET-2. We are deliberately 
agnostic about whether supermassive black holes or some 
other mechanisms provide the energy or momentum in these 
feedback prescriptions. For the first, we assume that major 
galaxy mergers ultimately drive the transformation of galax- 
ies from blue to red. We identify mergers in the simulation, 
and add a kinetic kick to the gas in the merger remnant to 
expel star-forming gas. This burst of ener gy mimics the ef- 
fects of a powerful quasar or starburst (e.g. lPi Matteo et al.l 
l2005h . 

In the second quenching model, we assume that the hot 
gas haloes surrounding massive galaxies facilitate quench- 
ing. We identify these hot, massive haloes in the simula- 
tion, then continually add thermal energy to the halo gas 
in order to cut off the fuel for star-formation. The addi- 
tion of th ermal energy mimics th e added heat from radio 
AGN (e.g. iMcNamara et al.ll2005 ) , unre solved gravitational 
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heating (e.g. 



Ijubelgas et al.l 



Pekel fc Birnboim |2008| ). conduction (e 
2004I ), or cosmic rays (e.g. iGuo fc Ohll200; 



2.1.1 Quenching via galaxy mergers 

Our first quenching prescription ties truncation of star for- 
mation to galaxy mergers. This requires identifying galaxy 
mergers within the simulations on the fly, then applying 
some quenching process to the particles in those mergers. 
We implement empirically-motived, high-velociy winds that 
eject cold gas from merger remnants. 

Merger identific ation is based on our on-the- fly galaxy 
finder, developed in lOppenheimer fc Pavel l|2008h . This al- 
gorithm uses a parallel friends-of-friends (FOF) algorithm 
to group together particles into galaxies. We set the link- 
ing length between particles to 0.06r a (H(z)/Ho) ll/3 , where 
r s is the mean interparticle separation, Ho is the hubble 
constant (70 km s -1 ), and H(z) is the hubbl e para meter as 
a function of redshift. lOppenheimer fc Pavel (2008) showed 
that stellar masses of groups found with this prescription 
successfully matched those obtained in post-processing with 
the slower but more physically-motivate d Spline Kernel In - 
terpolative Penmax (SKIP) Q algorithm l|Kere§ et al.ll2005l ). 
To optimize for speed, we divide the simulation volume into 
grid cells and only calculate distances between particles in 
neighboring grid cells. We identify groups independently for 
each central processing unit, then merge groups split onto 
adjacent CPUs by checking whether groups separated by less 
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Figure 1. Comparison of stellar mass histories calculated us- 
ing two different group-finders for three galaxies (one galaxy per 
panel). Dotted lines show normalised masses identified by SKID 
in post-processing, while solid lines show our faster on-thc-fly 
FOF masses offset by 0.75. Vertical line segments indicate merg- 
ers identified on-the-fly. (a) A galaxy with an easily identified 
major merger, (b) A galaxy with a major fly-by, incorrectly iden- 
tified as a merger, (c) A merger where the final coallescence is 
~ 1 Gyr later with SKID than with our on-the-fly method. The 
latter two problems are sufficiently rare to have little impact on 
our results. 



than 30 linking lengths indeed have particles that should be 
linked together. The frequency with which we run the galaxy 
finder varies from a few Myr to a few 10's of Myr, depending 
on the adaptive time-step intervals. 

We identify mergers on-the-fly by looking for persistent, 
rapid increases in stellar mass growth. This requires knowing 
the recent mass history of each galaxy. We track the 6 most 
recent galaxy catalogs constructed by the galaxy finder, then 
connect galaxies in the most recent catalog to their earlier 
most-massive progenitors. By tracking both the current and 
previous galaxy host of each particle, we can easily identify 
the progenitors of a given galaxy by searching through its 
constituent particles and noting their previous host galaxies. 

For a given galaxy, we identify a merger by compar- 
ing its current stellar mass to that of its recent progeni- 
tor galaxy. In detail, we compare the median mass of the 
three most recent progenitors (M reccllt ) to the median mass 
of the three older progenitors (M i dor ). If M reC ent /M u eI > 
(1 + 1/r) then we flag the galaxy as a merger. We use the me- 
dian approach, rather than just comparing the most recent 
galaxy mass to the previous one, to avoid declaring merg- 
ers of brief encounters and fly-by's. Since the usual value of 
r = 3 yielded good matches to the galaxy luminosity func- 
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tion in iGabor et al.1 |20ld ), we use the same choice here (see 
^14.41 for further discussion). 

This method of identifying mergers likely declares the 
merger before galaxies have fully merged. Since we use all 
particle types in the galaxy finder, two merging galaxies will 
appear as one once the central regions of their dark matter 
haloes overlap. Similarly, when two galaxies closely pass by 
each other, our method will declare a merger whether the 
galaxies eventually coalesce or not. Based on the work of 
iMaller et ail ((2006), this appears to be a generic problem 
with merger identification in cosmological simulations. We 
mitigate this effect partially by taking the median masses 
over several time steps as described above - rapid fly-by's 
will not alter the median mass enough to cause a merger 
identification, although some longer, close fly-by's do (as we 
show below). In cases where the two galaxies do finally co- 
alesce, the precise timing of the quenching does not matter 
much, since we effectively get rid of all cold gas in the rem- 
nant (see below). Early quenching could, however, stunt the 
star formation that might naturally occur during the merg- 
ing process, leading to slightly smaller stellar masses and 
lower metallicities in merger remnants. In the end, our re- 
sults are not expected to be sensitive to these effects. 

To illustrate the successes and challenges of our merger 
identification method, in Figure [T] we compare the stel- 
lar mass histories of galaxies identified on-the-fly to the 
masses of the same galaxies identified with SKID in post- 
processing. The mass growth histories from the on-the-fly 
FOF method broadly match those from SKID, and ~60% of 
galaxies with M ste ii a r > 10 10 Mq in this simulation undergo 
at least one merger in the last 10 Gyr. Panel (a) shows a well- 
matched single major merger. Panel (b) illustrates the diffi- 
culty with fly-by's that do not end up merging. Both FOF 
and SKID combine the two galaxies into a single galaxy at 
some point, and then separate them again as the two galax- 
ies move apart. Among galaxies with at least one merger, our 
method mis-identifies fly-by's as major mergers in < 10% of 
the cases. Panel (c) shows a case where detailed structure 
in the mass history identified by SKID is missed by FOF. 
In the FOF case, it looks like a clean-cut major merger. 
With SKID, however, the stellar mass jumps, then after a 
few hundred Myrs begins to fluctuate up and down before 
settling in a more massive state. This suggests an initial fly- 
by, then a few more close passages, then a final coallescence 
just over 1 Gyr after the initial merger identification. Among 
galaxies with at least one merger, our method identifies one 
of those mergers too early by > 500Myr in ~ 25%. In fu- 
ture work we may attempt to identify the final coallescence, 
and avoid mis-identifying fly-by's, by invoking a delay time 
of a few hundred Myr before actually declaring the merger. 
Given the broad nature of this study, such details of merger 
identification do not substantially affect our results. 

Once we have identified a merger, we apply ener- 
getic ejective feedback to the star-forming gas within the 
merger remnant, mimicking a powerful wind. Observations 
of post-starburst galaxies and QSOs suggest that such winds 
with velocities ~ 1000 km s _1 or more may be respon- 
sible for quenching star formation (|Tremonti et al.l 120071 : 
iFeruglio et afl|20ld ). We pursue the extreme case where all 
the star-forming gas is ejected in a wind because this is 
probably necessar y to produce a red galaxy remnant (cf. 
IGabor et al]|20ld ). Each gas particle in the merger remnant 
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Figure 2. Comparison of halo mass estimates from our on-the- 
fly algorithm (j/-axis) and a post-processing spherical overdensity 
algorithm (a-axis). The dotted line shows a perfect 1-to-l corre- 
spondence. In the on-thc-fly case, we multiply the baryonic (gas 
and stellar) mass found by our FOF galaxy finder by an empirical 
factor of 18 to get M es t with the best fit. This factor is larger than 
the expected cosmic factor (~ 6) because we tune the FOF link- 
ing length to match the stellar mass, not the total baryonic mass, 
of galaxies. Much of the baryonic mass remains in an extended 
halo not captured by the FOF method. The on-the-fly algorithm 
underestimates the halo masses of many satellite galaxies (blue 
crosses). This is because the spherical overdensity method asso- 
ciates satellite galaxies with their massive parent halo, whereas 
our FOF galaxy finder associates them with their smaller sub- 
haloes. 



galaxy is given a velocity kick of i>ki c k in the direction v x a, 
where v is the pre-kick velocity vector, and a is the acceler- 
atio n vector. 

iTremonti et~"aH l|2007t l compared post-starburst galax- 
ies with winds ^1000 km s" 1 to low-io nization broad ab- 
sorption line quasars l|Trump et al.|[200rJ) with winds up to 
~ 10 km s _1 . We choose Ukick = 1500 km s _1 as a com- 
promise for large-scale outflows that have likely exited the 
galaxy. Simulations with Ukick = 1000 km s _1 showed little 
difference in the results, since this is still well above the es- 
cape velocity for most haloes. As with our winds associated 
with star-formation, we decouple the wind particles from hy- 
drodynamic processes to enable escape from the ISM, just 
as we do with star formation-driven winds. 

As an addition to ejective feedback, we optionally im- 
plemented a thermal feedback mechanism associated with 
galaxy mergers. In this case, we heat all the star-forming 
gas particles to the virial temperature of the halo, as esti- 
mated below in t|2.1.2l It turns out that in combination with 
the ejective feedback, it makes little difference to the results, 
since ejected gas almost never returns to the galaxy whether 
we heat it or not. 



2.1.2 Quenching in hot haloes 

The presence of a hot gaseous halo around most massive red 
and dead ellipticals hints at some role in quenching. Hence 
for our second main quenching mechanism, we heat gas in 
the hot gaseous haloes around galaxies. For each galaxy 
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Figure 3. Hot gas fractions calculated on-the-fly vs. halo mass (M eB t) and spherical overdensity (SO) hot fraction for galaxies in a 
small volume simulation. Black dots show central galaxies and blue crosses satellite galaxies (as found using SO). Hot gas fractions are 
defined as the fraction of gas within the virial radius above 250,000 K. The left panel show s that low-mass h aloes have little hot gas in 
their surroundings, while increasingly larger haloes have increasingly high hot fractions (cf. iKeres et alj|2005h . The galaxies in low-mass 
haloes with f^ ot = 1 are satellite galaxies embedded in massive, hot haloes. The right panel compares our on-the-fly estimate to our 
post-processing, spherical overdensity estimate of /hot, and a thin solid line shows a perfect 1-to-l correspondence. We can reliably 
estimate the hot fractions of galaxies on-the-fly in the simulation. 



identified on-the-fly, we estimate the hot gas fraction in its 
halo. In haloes dominated by hot gas, we add heat by raising 
the temperature of gas particles around the galaxy. 

To obtain the hot gas fractions of galactic haloes, we es- 
timate the virial radius corresponding to each halo, find all 
the gas particles within that radius, then count up the mass 
on either side of a temperature threshold defining "hot." 
We estimate the halo virial radius, i? v ir, using the equation 
Pvir = Me 3 t/(4/37rJ?^ ir ). Here, p vlI is the virial density for 
collapsed structures at the appropriate epoch, and M cst is an 
estimated halo mass. We find Mest by scaling the baryonic 
mass of the galaxy from the group finder by a factor of 18. As 
shown in Figure [21 this empirical scale factor value best re- 
produces the total halo masses found using a more thorough 
spherical overdensity algorithm in post-processing. This de- 
parts significantly from the cosmic value of Om /fibaryon — 6 
because our on-the-fly group finder purposely does not cap- 
ture the full radial extent of the halo. Our FOF linking 
length is chosen to match galaxy stellar masses, so it is too 
small to capture much of the outer halo. Nevertheless, the 
empirical factor of 18 nicely estimates tht total halo mass. 

We add the mass of each gas particle within the virial 
radius to a tally of the hot gas mass if the particle's tem- 
perature exceeds 250,000 K; otherwise, we consider it "cold" 
or warm for the purpose of identifying hot gas haloes. We 
choose 250,000 (or 10 5 ' 4 ) K because this separates the main 
bulk of cooler IGM that eventually cools into galaxies from 
the hot tail of the IGM from which cooling and forming s tars 
is rare (cf. IKeres et al.1l2005l . l2009bl : iGabor et al]l20ld ). 

Once we have measured the hot and cold gas masses 
within the virial radius of each galaxy, we apply the quench- 



ing heat if / ho t = mhot/(m ho t + Jn co i d ) > 0.6. A hot frac- 
tion of 0.6 corresp onds roughly to halo m asses of 10 12 Mq 
i|Keres et al.ll2005l ). In IGabor etal] l|2010l ), we found that 
quenching star formation in haloes with of M > 10 12 Mq or 
/hot > 0.6 led to luminosity functions that best reproduce 
those observed. 

Figure [3] illustrates the effectiveness of our hot frac- 
tion estimate in a small simulation. Our on-the-fly hot frac- 
tion calcula t ion a ccurately reproduces the trend noted by 
IKeres et all l)2005l ): halos above ~ 1O 12 M have hot frac- 
tions above ~ 0.6 (left panel), and the transition between 
halos dominated by cold and hot gas is fairly abrupt. Out- 
liers from the trend with low masses and /hot = 1 are small 
satellite galaxies embedded in the hot halo of a larger galaxy. 
In the right panel, we compare our on-the-fly hot fraction 
estimate to that from our post-processing, spherical over- 
density method. With a high degree of reliability, we can 
determine whether the gaseous halo surrounding a galaxy is 
hot. 

In FOF galaxies above the hot fraction threshold, we 
raise gas particles to the estimated virial temperature of the 
halo. In our base model, we only raise the temperatures of 
particles that are not star-forming, i.e. that are not within 
the existing cold ISM of galaxies within the halo. We also 
do not raise the temperature of particles alread y exceeding 
the virial temperature. We use equation 59 from lvoitl (|2005r ) 
to estimate the virial temperature, using M es t for the halo 
virial mass: 

Here, ks is Boltzmann's constant, H(z) is the redshift- 
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dependent Hubble parameter, and Ho is the Hubble con- 
stant (at z = 0). In practice, we found that instantaneously 
heating particles to the virial temperature sometimes led 
to unphysically high hydrodynamic accelerations. To pre- 
vent adaptive time-steps from becoming too short, we in- 
stead double the entropy of the gas particle until it reaches 
the virial temperature. Through this heating mechanism, we 
prevent gas from ever cooling and forming stars. 

A key aspect of our mechanism is that energy is contin- 
ually added to keep particles hot. Gas within the hot halo 
that cools below T v ir, or accreted gas that is below T v i r , is 
heated as described above during every time step of the sim- 
ulation. This is perhaps overkill; to suppress star formation, 
it may not be necessary to keep gas at or above the virial 
temperature, nor to do so at every timestep. Here, however, 
we will make this maximal heating assumption to ensure 
star formation is suppressed, even though we will discuss 
later that the energetics of this are quite extreme. We will 
explore relaxing the heating assumption in future work. 

Given the overall success of our hot gas quenching 
model, we explored many variants on our base model: 

(i) We heat all gas within halos to the virial tamperature, 
including star-forming gas within the ISM of galaxies. 

(ii) We use a criterion based directly on the halo mass 
(rather than hot gas fraction). While there is a close cor- 
relation between hot gas haloes and massive dark matter 
haloe s (cf. Figure H iBirnboim fc Dekell 120031 ; | Keres et al.l 
2005), some models have had success quenching star forma- 
tion above a particular ha l o mass, typically a round 10 1 Mq 
(e.g. ICattaneo et alJl20o5lGabor et alj|20ich . 

(iii) We restrict the heating to within a certain radius 
^heat from the center of the halo. Thi s is an attempt to mimic 
feedback from a central radio AGN. IVoit fc Donahue! (|2005h 
argue, based on observations of an entropy f loor within the 
centra l 30 kpc in cooling flow clusters from iDonahue et al.l 
(2006), that an AGN outburst sufficient to counteract the 
cooling will create subsonic, buoyant bubbles at a typical 
radius of 30 kpc. The region outside this radius should be 
heated by the expanding bubbles. Following this, we use 
^heat = 30 kpc. A numerical difficulty with this is that 
we found that large galaxies in our simulations have star- 
forming regions at significantly larger radii. This result oc- 
curs for galaxies identified using both FOF and SKID meth- 
ods. This may owe either to satellite galaxies at these radii 
being grouped into the central galaxy, or to the threshold 
density for star formation being too low. We note that the 
SF threshold is « 0.1 cm -3 , which is well below ISM densi- 
ties where star formation occurs in reality. This is a histori- 
cal numerical convenience motivated by the inability of our 
simulations to fully resolve ISM densities. In any case, this 
variant produces a red sequence, but has too many massive 
galaxies remain blue due to the ongoing star-formation. Ap- 
parently in our simulations, heating only the central regions 
of a halo does not heat the outer regi ons as suggested b y 
observations of sound waves in clusters (Fa bian et al.| [2006'). 
The failure of central heating to spread energy to large radii 
may suggest that we do not adequately resolve the hydrody- 
namics in our typical hot haloes, or that raising the temper- 
ature to only Tvir may be insufficient, and we need to raise 
it further. 

(iv) We do not heat gas below our 250000 K critical tem- 



perature. This is intended to allow for cold flows from IGM 
fil aments to penetrate into hot haloes, as has been argued 
by iDekel et ail l|2009bl ). We discuss this variant further in 

am 

(v) We do not heat gas within sub-haloes. As indicated 
by Figure [21 many small satellite galaxies live in the hot 
gas haloes of their larger hosts. AGN in these small galaxies 
may or may not contribute significantly to the overall heat- 
ing of the hot gas halo. Since their importance in quenching 
massive galaxies is debatable, we implemented a method to 
identify these sub-haloes and exclude them from the feed- 
back prescription. When calculating hot gas fractions as de- 
scribed above, we estimate the virial radius of each FOF 
galaxy. If a galaxy's center falls within the virial radius of 
another, more massive, galaxy, then we consider it a sub- 
halo. Once we identify sub-haloes, we do not apply the ther- 
mal feedback to gas within them. We discuss results of this 
mechanism in N4.3I 

2.2 Simulation parameters 

Simulating massive galaxies presents a substantially larger 
computational challenge than simulating more common sys- 
tems. One is pushed towards larger volumes (> 50 Mpc/h) 
by the need to obtain a sufficient statistical sample of mas- 
sive galaxies. Conversely, since massive galaxies begin form- 
ing stars first, this requires simulations that robustly resolve 
early star formation, requiring higher resolution (~few kpc). 
Given computational resource limitations, we can only per- 
form a few large runs. To explore all the variants listed 
above, we are forced to run smaller simulations. We choose 
to compromise on volume to ensure robust star formation 
histories. 

Our primary simulations use a box size of 48ft _1 Mpc 
and a gravitational softening length of 3.75h~ 1 kpc (equiv- 
alent Plummer), with 256 3 dark matter and 256 3 gas parti- 
cles. This yields a gas particle mass of 1.2 x IO 8 M0, and a 
dark matter particle mass 5.1 x larger. Our smaller simula- 
tions keep the same mass and spatial resolution, but reduce 
the box size to 24/i _1 Mpc with I28 3 dark matter and I28 3 
gas particles. All simulations and analyses use a Wilkin- 
son Microwave Anis otropy Probe concordance cosmology 
(|Komatsu et alj|2009h with Ho = 100ft = 70km s" 1 Mpc" 1 , 
matter density f2 ro = 0.28, baryon density Qt = 0.046, a 
cosmological constant with Qa = 0.72, root mean square 
mass fluctuation at separations of 8 Mpc as = 0.82, and a 
spectral index of n — 0.96. Our small runs are useful for 
testing some of the variations in the quenching physics, and 
also for tracking details of each on-the-fly galaxy and feed- 
back event. The larger simulations do a much better job of 
sampling the massive end of the luminosity function, but 
take many weeks to run. 

We focus primarily on three models which we refer to 
as no quenching, merger quenching, and hot gas quench- 
ing. The no quenching model includes star formation and 
momentum-conserving winds, but no additional explicit 
feedback mecha nism. It is effectively th e same model as the 
"vzw" model of lOppenheimer fc Pavel l|2008T ). The merger 
quenching model incorporates 1500 km s winds from 
merger remnants to quench star formation. For the hot gas 
quenching, we focus on a model where we heat only the non- 
star-forming gas particles in FOF galaxies with /hot > 0.6. 
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In f|4]we will also explore some simulations with slight vari- 
ations of the quenching models. 



2.3 Simulation outputs and analysis tools 

We output snapshots of each simulation at 108 redshifts, 
starting at z — 30 and ending at z = 0. The snapshots 
contain information for every simulation particle, including 
position, velocity, mass, metallicity, gas density, gas temper- 
ature, star-formation rate, and time of formation (for star 
particles). In addition to these outputs, we save the galaxy 
catalog from each instance of the on-the-fiy galaxy finder for 
later analysis. 

Our suite of analysis tools allows us to compare our 
simulation results directly with observations. SKID provides 
a list of member particles (star and star-forming gas) for 
each simulated galaxy. The sum of star particle masses is 
then the galaxy stellar mass, and the star formation rates 
of the gas particles are summed to give the star formation 
rate of the galaxy. 

We then calcu l ate ga laxy spectra using the models of 
iBruzual fc Chariot! (|2003l ). We treat each star particle as a 
single stellar population with an age and metallicity deter- 
mined directly in the simulation. By adding up the spectra 
of all star particles within a galaxy, we obtain the spec- 
trum of that galaxy. Then we convolve galaxy spectra with 
observational filter curves to obtain broad-band colors and 
magnitudes. 
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3 OBSERVATIONAL CONSTRAINTS 

To characterize the effectiveness of our quenching models, 
we compare the results of our simulated galaxy populations 
to galaxy observations in the low-redshift (z < 0.1) uni- 
verse. Our goal is to create a realistic red sequence, so we 
focus on color-magnitude diagrams (CMDs) and red-galaxy 

luminosity funct ions (LFs). 

We follow iGabor et all l|2O10D in using the low- 
redshift version of t he Value-Add ed Galaxy Catalog (VAGC, 
Blanton et"ai]|2005h of the SPSS (lAdelman- McCarthy et all 



20081 ; iPadmanabhan et all I2008T ). This 



catalog includes 
ugriz + JHK absolute magnitudes for ~ 170000 galaxies 
with z < 0.05. After converting absolute magnitudes to our 
preferred cosmology (with h = 0.70), we can straightfor- 
wardly plot CMDs. We will plot g— r vs r because the r-band 
light is a good tracer of stellar mass. For the luminosity func- 
tions, we use r-band luminosities and the 1/V max method 
ijSchmidtll 19681 ) with the Vmax values in the VAGC. We com- 
pare these observational results with simulation snapshots at 
z — 0.025, which w e will refer to as low redshift or z — 0. 

As discussed in IGabor et all (|20ld ) and £|4.2l uncertain- 
ties in dust obscuration motivate us to compare mass func- 
tions (MFs) of blue galaxies rather than luminosity func- 
tions. For this we use publicly available stellar masses de- 
rived for SPS S galaxies by fitting obs erved galaxy spectra 
to templates (jKauffmann et al1l2003bl ). We cross correlate 
this sample with the VAGC to obtain ~ 40000 galaxies with 
stellar masses and broad-band colors for the MFs. 



Figure 4. CMD (top panel) and LF (bottom) for a simulation 
with no quenching. In the CMD, grayscale represents the density 
of SDSS galaxies, points are simulated galaxies, the dot-dashed 
line is a fit to the observed red sequence, and the solid line sepa- 
rates the simulated red sequence and blue cloud. Without quench- 
ing, our simulations do not produce red galaxies. In the LF, sym- 
bols represent the observed luminosity functions for red and blue 
galaxies, and lines are for simulated galaxies. The red galaxy lu- 
minosity function is too low by roughly an order of magnitude. 
Simulated blue galaxy luminosities are computed without a cor- 
rection for dust to distinguish the intrinsic red sequence more 
easily. 



4 CONSTRAINING QUENCHING MODELS 
4.1 Basic models: Red Galaxies 

Since our goal is to reproduce a realistic red sequence, we 
compare the red galaxy CMPs and LFs produced by our var- 
ious models to determine the best match. We first illustrate 
the problem with simulations that do not have additional 
quenching mechanisms. Figure [4] shows the full CMP and 
LF for a simulation with no quenching feedba ck, compared 
to observations. Following IGabor et al] l|2010h , we separate 
red simulated galaxies from blue in the CMP using the solid 
line. This line is bluer and shallower than that used for ob- 
served galaxies because it better captures the bimodality 
in our most successfu l models. As we discuss below (and 
m IGabor et ai]|201Ch , our model red sequence is generally 
too blue and has little or no slope. We estimate LF uncer- 
tainties based on a jackknife re-sampling of the simulation 
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Figure 5. CMDs zoomed in on the red sequence (top 3 pan- 
els) and LFs (bottom panel) for simulations with no quenching, 
merger quenching, and hot gas quenching. This shows the main 
result of this work: merger quenching fails to produce a substan- 
tial red sequence, whereas hot gas quenching succeeds. 



volume |Finlator et al.ll2006l ) (typical observational uncer- 
tainties are smaller than the symbols except for the brightest 
bins). 

Almost all galaxies in this no-quenching simulation con- 
stantly accrete gas from the IGM, and this gas supplies 
the fuel for ongoing star formation. As a result, the sim- 
ulation produces few red galaxies and many very massive, 
very bright blue galaxies. The LF reflects this effect. The red 
galaxy LF is too low (relative to that observed) by an or- 
der of magnitude for r ~ —21.5, and the blue LF is too 
high. Furthermore, this model does not produce the ob- 
served sharp turnover in the lum inosity function at high 
masses (|Qppenheimer et"al]|2010l ). 

We show the simulated blue galaxies without account- 
ing for obscuration from dust to illustrate the intrinsic 
galaxy luminosities. Dust obscuration will push the blue 
galaxies up and to the right (redder and fainter) in the CMD, 
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Figure 6. Cosmic histories of mass (solid lines) and instanta- 
neous star-formation rate (dotted lines) for 4 galaxies selected 
from a simulation with merger quenching. These illustrate the 
variety of outcomes following major mergers. Vertical solid lines 
above mass histories indicate major mergers, and we show the 
2 = stellar mass (MJ in Mq ) and maximum star- formation rate 
(in MQvr - 1 ) for each, (a) A galaxy with moderate, fairly con- 
stant star-formation rate, (b) A galaxy with some major mergers. 
By design, the star-formation rate drops to zero immediately fol- 
lowing a merger. The star-formation rate recovers within about 
2 Gyrs due to gas accretion and possible minor mergers, (c) A 
galaxy whose major merger quenches almost all subsequent star 
formation until z = 0. (d) The complex mass history of the central 
galaxy in one of the two most massive groups in the simulation. 
This galaxy undergoes several major mergers, after each of which 
the SFR recovers within 1 Gyr due to fresh inflows of gaseous 
fuel. 
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Figure 7. Star formation histories (dotted lines) and stellar 
masses (solid lines) for three galaxies in our hot gas quenching 
model (cf. Figure |6j. (a) A galaxy with fairly quiescent star- 
formation and steady mass growth. At z S3 0.3, this galaxy's halo 
becomes dominated by hot gas, after which we begin heating the 
gas to prevent cooling. Starved of new fuel, star-formation de- 
clines to zero over ~ 2 Gyr. (b) A minor merger (lookback time 
~ 2 Gyr) delivers fuel to a galaxy which had begun quenching 
about 2 Gyr earlier, (c) A major merger increases a galaxy's halo 
mass to the point where it can support a hot halo, after which it 
is quenched. 



possibly contaminating the red sequence. The observed uni- 
verse does not contain star-forming galaxies as massive as 
those shown here, so estimating the extinction in these sim- 
ulated galaxies would require large extrapolations that are 
unlikely to be robust. We discuss dust extinction and blue 
galaxies further in i|4.2l 

Focusing just on the red sequence, we compare results 
for models with no quenching, merger quenching, and hot 
gas quenching in Figure [5] Like the no quenching simula- 
tion, the merger quenching model leaves a very sparse red 
sequence. Our hot gas quenching model, however, produces a 
red sequence of galaxies well-separated from the blue cloud. 

The luminosity function panel in Figure [5] represents 



one of the central results of this paper. It emphasizes that 
our hot gas quenching model nicely produces a sharp trunca- 
tion in the red galaxy LF, while our merger quenching model 
does not. Because it produces so few red galaxies, merger 
quenching does not significantly improve the red galaxy lu- 
minosity function of our simulations without quenching. The 
hot gas quenching model fares well, matching the observed 
red sequence luminosity function at all magnitudes we re- 
solve. 

For the very brightest galaxies, the relatively small vol- 
ume of our simulations hampers the comparison with data - 
we do not effectively sample the highest peaks in the density 
distribution that lead to the most massive galaxies. This is 
reflected in our errorbars. 

Examining the halo quenching red sequence further, we 
see that it is too blue by about 0.1 magnitudes in g — r, and 
has a nearly flat slope that fails to match the observed trend 
that brighter red galaxies a re red der. We examined this issue 
at length in iGabor et all {2010), and the same underlying 
causes are likely here. The overall "blueness problem" may 
result from either metallicities that are too low or ages that 
are too young. Our red galaxies have metallicities lower than 
those inferred for the observed red sequence in the VAGC, 
and lower metallicities lead to bluer stars. Mean stellar ages 
of our galaxies are at least as old as those inferred from 
the observations, so incorrect ages cannot cause the red se- 
quence to be too blue. Since we match the stellar masses, the 
metallicity deficit is unlikely to owe to an underproduction 
of stars. It could be that our outflow model ejects too many 
metals from these systems, or else we might have adopted 
supernova metal yields that are too low. The red sequence 
slope poses a subtler problem, perhaps suggesting that sim- 
ulated massive galaxies merge with too many, too massive, 
or too young satellite galaxies with lower metallicities. In 
i]4.2l we describe a metallicity re-calibration that matches 
our simulated red galaxy colors to those observed. 

Why does merger quenching fail? In lGabor et al.l (|2010l ) 
we showed that if star formation is completely quenched in 
the remnants of major mergers, then our models can produce 
a reasonable match to the red sequence. Merger quench- 
ing therefore does not fail because there are too few major 
mergers. Rather, merger quenching fails because star for- 
mation resumes even after a merger remnant loses all its 
star-forming gas. 

We illustrate the resumption of star formation after 
merging in Figure [U Each panel corresponds to an individ- 
ual galaxy taken from a small test simulation with merger 
quenching. For each galaxy we show the normalized star- 
formation history (dotted lines) and stellar mass history 
(solid lines) over cosmic time, with vertical tick marks in- 
dicating major mergers. Galaxies with no major mergers 
(panel a) have fairly constant star formation histories (since 
z ~ 2), with peaks and valleys depending on the details 
of accretion from the IGM and outflows powered by stars. 
Galaxies that undergo major mergers (panels b-d) generi- 
cally stop forming stars when a merger occurs because we 
instantaneously eject all the cold gas. However, in every case 
the merger remnant resumes star formation as new fuel is 
accreted from the IGM. The timing and intensity of the re- 
sumption of star-formation varies, probably depending on 
environment. A very massive galaxy living at the center of 
a large potential well (panel d) will accrete new gas and re- 
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turn to high star-formation rates in ~ 1 Gyr. Some smaller 
galaxies (e.g. panel c) never return to SFRs near their pre- 
merger levels, though they do resume some star-formation 
within ~ 2 Gyrs. 

As a comparison, we show star-formation histories for 
quenched galaxies in our hot gas quenching model in Figure 
[7] The galaxy in panel (a) grows steadily in stellar mass 
until its halo becomes dominated by hot gas (at z « 0.3). 
At this point our hot gas heating cuts off additional fuel 
for star-formation, and the galaxy gradually uses up its pre- 
existing reservoir of gas. Star formation peters out after ~ 
2 Gyrs. Even after the starvation process begins, merging 
galaxies can momentarily supply new fuel, as reflected by a 
jump in SFR in panel (b) at a lookback time of ~ 2 Gyrs. 
In panel (c) we show a galaxy whose major merger raises 
its halo mass enough that a hot gas halo becomes stable. 
These star-formation histories all show a slow decline in star- 
formation rate over ~ 2 Gyr timescales, similar to the star- 
formation timescale in our sub-grid star-formation model 
|Springel fc Hernquistl20"03l ). After star formation stops, the 
galaxy stellar mass declines slowly owing to stellar evolution. 

In summary, a simple model where all star formation is 
quenched in galaxies with a hot halo fraction above 60% pro- 
duces a red galaxy luminosity function in very good agree- 
ment with observations. It also produces a reasonable red 
sequence, althou gh it is still too blue and too shallow as we 
saw previously in lGabor et al.l (|2010l ) , likely reflecting issues 
with enrichment. In contrast, assuming that all major merg- 
ers are associate d with 1500 km s" 1 outflows as suggested 
by observations l|Tremonti et al.ll2007T ) does not significantly 
populate the red sequence. 



4.2 Basic models: blue galaxies 

We have shown that additional heating in hot gas halos can 
produce a red sequence whose luminosity function matches 
that observed in the local universe. We now examine galaxy 
colours more closely. Since these are sensitive to uncertain 
dust corrections and problems with metallicities, we will use 
galaxy stellar mass functions later in this section. 

Figure [8] shows a complete CMD and corresponding LFs 
for our hot gas quenching model. Here, we apply an ad hoc 
metallicity correction to force the simulated red sequence to 
match the observed colors. This assumes that our predicted 
stellar ages are correct, but that our simulations underpre- 
dict stellar metallicities in massive red galaxies. We add a 
metallicity of Z a( jd to each single stellar population (i.e. star 
particle) when calculating galaxy colors using stellar popu- 
lation synthesis models. By trial and error, we find a good 
fit using the relation 



0.003+ (6 x lO- 14 M stc n ar /M ). 



(2) 



These are absolute met al mass fractions, an d solar metal- 
licity here is « 0.012 l|Asplund et all 120051 ). This correc- 
tion gives an indication of how far off our metallicities must 
be. For a 10 11 Mq galaxy (with typical absolute magnitude 
r ~ —22), Z a dd = 0.009, representing an increase of ~ 50% 
over typical simulated metallicities for such galaxies. 

We show blue galaxies in the CMD with a correction 
for dust extinction of optical light in the g and r bands. 
In our dust prescription, the E(B — V) extinction scales 
with star-formation rate based on empirical relations with 
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Figure 8. Full CMD (top panel), LFs (middle panel), and stellar 
mass functions (lower panel) for our hot gas quenching model. The 
simulated CMD includes a metallicity re-scaling to best reproduce 
the observed red sequence colors, as well as dust reddening and 
extinction for the blue galaxies. Due to uncertainties in the dust 
prescription, we prefer to use stellar mass functions (lower panel) 
to compare simulated and observed blue galaxies. Our hot gas 
quenching model yields too few massive blue galaxies. 



UV luminosity of loc a l galaxies dWang fc Heckmanl 1 19961 ; 
I Somerville et all l20t)ll ; iFinlator et al.l 120061 ). We assume a 
ICalzetti et al.l (2000) reddening law. This dust prescription 
yields a blue cloud in fairly good agreement with data, espe- 
cially in the qualitative shape. It also shifts a small number 
of intrinsically blue galaxies onto the red sequence, but we 
count these galaxies as blue for the purposes of calculating 
the lu minosity and mass fu nctions. 

In iGabor etafl (|2010h . we discussed the difficulties in 
modeling extinction. Given a model galaxy with a speci- 
fied star-formation and gas accretion history, we need to 
know the level of dust extinction and reddening to com- 
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Figure 9. Stellar mass functions split into blue and red galax- 
ies (colors) for a simulation with no quenching (dashed lines), a 
simulation with merger quenching (solid lines), and SDSS data 
(triangles). Inducing powerful outflows after galaxy mergers does 
suppress the growth of massive blue galaxies, but it does little to 
increase the number of red sequence galaxies above that produced 
in the no-quenching simulation. 



pare colors with observed galaxies. Several dust extinc- 
tion prescriptions keyed to physical properties of galax- 
ies, such as metallicity or sta r-formation rate, have been 
proposed l|Finlator et al.ll2006l ). Different dust prescriptions 
yield substantially different colors and magnitudes for our 
star-forming model galaxies, suggesting large uncertainties. 
This is especially true for large, bright blue galaxies since 
these are not well represented in the real Universe. For red 
galaxies th is is not a problem because they have very little 
extinction l|Lauer et al.ll2005h . 

One way to mitigate the large uncertainties in dust 
obscuration is by examining stellar mass functions for the 
blue galaxies, rather tha n luminosity funct i ons. W e use stel- 
lar mass estimates from lKauffmann et alj (|2003bl) based on 
SDSS spectra, which account for extinction self-consistently 
within the SED fitting procedure for each galaxy. We cross- 
correlate these stellar mass measurements with galaxies from 
the VAGC, allowing us to use broad-band colors and gen- 
erate mass functions using the l/Vmax method. Using the 
broad-band colors, we separate the galaxies into the red se- 
quence and blue cloud as before. 

We show the resulting stellar mass functions for our hot 
gas quenching model in the bottom panel of Figure [8] As ex- 
pected from the luminosity function results, and given the 
tight correspondence between luminosity and mass for old 
stellar populations, the red galaxy mass function matches 
the observations well. The blue galaxy stellar mass function 
is overall a good match to data as well, although it shows 
a slight underproduction of star-forming galaxies just above 
the knee of the LF. The small number of galaxies (< 10) 
populating the bright end show only trace levels of star- 
formation. Our hot fraction criterion (/hot > 60%) for feed- 
back results in a fairly sharp critical halo mass of - 1O 12 M 
(Figure [3]), and therefore stellar mass of ~ 1O 1O ' 5 M , above 
which fueling of star-formation is turned off. Nevertheless, 
later infall of blue systems can turn a massive galaxy some- 



what bluer temporarily. Hence this model produces massive 
blue galaxies, although the details of that are quite sensi- 
tive to our quenching prescription and should be taken with 
some caution. 

This quenching model minimally impacts small galax- 
ies, basically on l y whe n they are satellites within larger ha- 
los. iDave et all l|201lf ) show that the fraction of satellite 
galaxies at these masses is about one-third, which explains 
why any impact on satellites results in only minor changes 
to the mass function. At the faint end, our model some- 
what overproduces the faint star-forming galaxy population, 
which i s a result also seen in the total mass function exam- 
ined in lOppenheimer et al.l (|201dh ■ Thus, while our hot gas 
quenching model successfully matches the observed number 
densities of red sequence galaxies, it fails to match all the 
details of the blue galaxy population. 

For comparison, we show the galaxy stellar mass func- 
tions for our no-quenching and merger outflows simulations 
in Figure [9] Solid lines correspond to our no-quenching 
model, and dashed lines to merger quenching. By ejecting 
the gas from galaxy merger remnants, we substantially re- 
duce the number density of very massive star-forming galax- 
ies. As indicated by the luminosity function analysis above, 
merger quenching does not substantially increase the num- 
ber of red galaxies. In both models, both red and blue mass 
functions look roughly like power laws, lacking the expo- 
nential cutoff observed in SDSS. This shows that merger 
quenching helps with only one of the two problems in the 
simulated LFs, namely that it reduces the number of blue 
galaxies, but it does not increase the number of red galaxies. 

In summary, quenching by halo hot fraction successully 
reproduces the SDSS red galaxy LF, and broadly repro- 
duces the SDSS red sequence. An empirical augmentation 
to galaxy metallicities given in equation [2] which increases 
metal content by typically 50% (with a range of 10-100%), 
results in an excellent match to the red sequence. The pre- 
dicted red galaxy mass function is likewise an excellent fit to 
data, and the blue galaxy mass function broadly matches as 
well although there is a hint of too few bright galaxies and 
too many faint ones. The agreement using hot halo quench- 
ing is far superior to that using merger quenching, as the lat- 
ter does not increase the red sequence over a no-quenching 
model although it does significantly suppress blue galaxies. 
These simulations therefore suggest that halo gas quenching, 
by itself, drives the formation of the red sequence. 

4.3 Model variations 

We now consider physically-motivated variations on our ba- 
sic models. Here we present three variations that all include 
hot gas heating. In the first we combine our merger winds 
quenching model with our hot gas heating model. Both of 
these mechanisms are well-motivated by observations, hence 
both processes may occur in the real universe. We refer to 
this as the "combined" model. 

In a second variation, we add an exception to the hot 
gas heating model. Several authors have argued that cold 
flows of gas along filaments can penetrate hot gas haloes at 
high redshift, since filaments are denser at earlier epochs. 
In principle our simulations should capture this behaviour 
with their hydrodynamic treatment. Our basic hot gas heat- 
ing model, however, bypasses the hydrodynamics to heat all 
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Figure 10. CMDs (top 3 panels) and LFs (bottom panel) for 
three alternative quenching models: combined hot gas heating + 
merger winds (top CMD, dotted line in LF), the same combined 
model except without heating any gas below 250000 K (middle 
CMD, dashed line in LF), and a model with hot gas quenching 
without heating gas around sub-haloes (bottom CMD, solid line 
in LF). The combined model performs well, with slightly lower 
LFs than the hot gas quenching-only model. Allowing cold flows 
fails due to wind recycling, and allowing gas fueling in sub-haloes 
leads to discrepancies at both low and high luminosities in this 



(non-star-forming) gas around a galaxy. In order to prevent 
artificial heating of these cold flows, allowing them to pen- 
etrate hot haloes, we implemented a model where we only 
heat gas that is already above our hot gas cutoff of 250000 
K. We present this variation with the combined mergers + 
hot gas model, but we note that the results are similar when 
we apply this variation to our basic hot gas heating model. 
We refer to this as the "cold flows" model. 

In a third variation, we tested a model where we do not 
heat the hot gas in the vicinity of sub-haloes. In our simu- 
lations, some sub-haloes are embedded in their parent's hot 
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Figure 11. Blue and red galaxy mass functions (color-coded) for 
our combined quenching model (solid lines) and our model with- 
out heating gas around sub-haloes (dashed lines). The combined 
model shows a larger deficit of massive blue galaxies, and the sub- 
haloes model shows a similar deficit to that in our basic hot gas 
heating model. 



gas halo, even if the sub-halo is too small to support its own 
hot gas. Our basic quenching model heats the gas in these 
sub-haloes just the same as in the massive, parent haloes. 
Heating from AGN feedback is typically observed in the cen- 
ters of groups and clusters, and it is not clear whether heat- 
ing from AGN in small satellite galaxies actually occurs in 
the real Universe. In our simulations, many satellites partic- 
ularly at low mass es are stripped o f gas owing to ram pres- 
sure or starvation (|Dave et al.ll201ll ) , but larger satellites are 
typically star-forming. X-ray observations do suggest that 
most satellites in clusters have mini cooling cores analogous 
to those seen in brightest cluster galaxies, and that heating 
by supernovae o r AGN may be necessary to quench them 
l|Sun et al.ll2007h . Some s atellite elliptical galaxies show evi- 
dence for AGN outbursts (jjones et al.| [2002 ; Mac hacek et al.l 
l200fj) . 

To distinguish sub-haloes (or satellite galaxies) from 
their parent haloes, we use the virial radii estimated when 
calculating hot gas fractions. If we consider two haloes A 
and B, then halo A is a subhalo of halo B if and only if (1) 
the center of halo A falls within halo B's virial radius and (2) 
halo B is more massive than halo A. We flag sub-haloes in 
our group catalog, and we do not add heat to their gas. Note 
that in this "subhaloes" variant, we do not induce outflows 
in galaxy merger remnants. 

We show results from these three model variations in 
Figures [101 and [TT1 Figure [TOl analogous to Figure [5] shows 
CMDs and LFs for just the red sequence in our three variant 
quenching models. Combined quenching has only a small 
overall effect on the red sequence when compared to our 
basic hot gas quenching: at all masses, the number density of 
red galaxies is slightly below that observed. In short, merger 
quenching has little overall effect on the galaxy population. 

Our "cold flows" model fails to produce a red sequence. 
This is not due primarily to pristine infall along filaments, 
but rather due to recyc ling of material expelled by star 
formation-driven winds (|Oppenheimer et all |2010| ). Wind 
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recycling dominates infall in galaxies since z ~ 1, and is 
stronger in more massive systems that have short recycling 
times. Although particles in these systems are ejected at 
speeds of up to ~ 1000 km/s, their cooling times times are 
short owing to the high level of enrichment, and so many 
never heat above the shock threshold. We discuss this issue 
further in ij6] 

Our variation that excludes heating in sub-haloes (the 
"subhaloes" model) increases the number of bright red 
galaxies and decreases the number of faint (r < —20) 
red ones. The deficit of faint red galaxies is straightfor- 
ward - turning off quenching in satellites allows those satel- 
lites to accrete gas and remain blue. The excess of bright 
red galaxies arises because the massive galaxies that make 
up the bright end grow significant mass via mergers with 
their satell ites, even af t er th e central galaxies have been 
quenched l|Gabor et al.l [2010). Without quenching, these 
merging satellites form more stars before they are subsumed 
than they do with quenching. These mergers then result in 
more massive central galaxies, and a corresponding increase 
in the bright-end luminosity function. 

To study the blue galaxies in these variant models, we 
again look at stellar mass functions. Figure [TT1 shows mass 
functions for the combined quenching model (solid) and the 
hot gas quenching model without sub-halo heating. In com- 
parison to our basic hot gas heating model, the combined 
model results in a sharper cut-off in the blue galaxy mass 
function above M stc ]iar ~ 10 10 ' 5 Mq, leading to a larger 
deficit of massive blue galaxies. The model without sub-halo 
quenching yields a similar blue galaxy mass function to our 
basic hot gas heating model. 

Overall, the variant models we tested do not solve the 
underlying problems with our basic models, while introduc- 
ing new problems. Our combined model exacerbates the 
deficit of blue galaxies at masses > 10 10 ' 5 M . Our model 
with cold flows fails to produce a red sequence. Our model 
without heating of sub-haloes creates discrepancies in the 
bright and faint end red galaxy LF. 

The combined model leads to similar result s as the basic 
model , but with fewer red sequence galaxies. In lGabor et al.l 
|2010l ) we showed that the halo mass cutoff, or correspond- 
ingly the hot fraction cutoff, above which we quench star- 
formation largely determines where the exponential cutoff 
begins in the LF. We therefore expect that raising fj££ from 
0.6 to, say, 0.7 might bring the combined quenching model 
into better agreement with observations. We cannot there- 
fore easily distinguish between the basic hot gas quenching 
and the combined quenching models. 

In summary, by testing many variants on our basic 
quenching models, we find that hot gas quenching is re- 
quired to get a red sequence, while merger quenching has 
minimal impact. Variations based on not quenching gas in 
cold flows or subhalos have a noticeable but not large ef- 
fect on the LF and CMD. None of these variants impact the 
slope or blueness of the red sequence, which continue to be 
a difficulty. Hence while all these quenching model variants 
are physically motivated, in the end a model with simple 
hot gas quenching matches the ensemble data as well as or 
better than any of the others. 



4.4 The dependence of merger quenching on the 
critical mass ratio 

Studies typically define major mergers as those with a 3:1 
mass ratio or smaller. We have adopted this critical mass ra- 
tio (r = 3) here for merger quenching, assuming that merg- 
ers with larger mass ratios have no effect. In this appendix, 
we show that allowing larger mass ratio mergers to trigger 
quenching does not lead to a red sequence in better agree- 
ment with observations. Although using r > 3 leads, in a 
general sense, to more quenching, it does not lead to a bi- 
modality in galaxy colors nor a significant number of red 
galaxies. 

Figure [12] compares color-magnitude diagrams and lu- 
minosity functions for four different simulations: one with 
our preferred hot gas quenching model, and our merger 
quenching model with three values of the critical mass ra- 
tio, r — 3, 5, and 10. These simulations are run with the 
same resolution as our main simulations, but with a vol- 
ume 8x smaller. While this volume does not sample mas- 
sive galaxies well, it does illustrate the main trends. In par- 
ticular, we show the hot gas model because it exhibits a 
well-populated red sequence distinct from the blue cloud. 
We show marginally resolved galaxies with r— band > — 20 
to help clarify the underlying behavior of the simulations. 

As described in the text, the merger quenching model 
under-produces red galaxies while suppressing overall galaxy 
growth (compared to a model without quenching). This is 
true for any value of the critical mass ratio. Furthermore, the 
sparsely populated region of red galaxies is poorly separated 
from the blue cloud, representing the tail of a unimodal dis- 
tribution rather than one peak of a bimodal distribution. As 
we move to higher values of the mass ratio, overall galaxy 
growth is suppressed more and more, but only a few galax- 
ies at any given time are red. Even with frequent quenching 
events due to minor mergers, galaxies continuously accrete 
new material from the IGM to fuel star-formation. 



5 PHYSICAL IMPACT OF QUENCHING 

Our quenching models have many physical and observable 
implications for galaxies and the IGM. Winds from merger 
quenching could impact the diffuse IGM. Our hot halo 
quenching has a major impact because it changes star for- 
mation histories. Since outflows associated with star forma- 
tion enrich the IGM and often end up returning to galaxies, 
changes in the star-formation rates can lead to changes in 
IGM and galaxy properties in complicated ways. Exploring 
all these changes is beyond the scope of this work, but here 
we explore two physical consequences of quenching - feed- 
back energetics, and the overall build-up of the red sequence. 



5.1 Quenching energetics 

In our simple implementations of quenching mechanisms, we 
have made no effort to limit the energy available to, e.g., that 
accreted onto a supermassive black hole. We simply eject or 
heat up the prescribed gas particles. Here we investigate 
whether the energy added via quenching is reasonable and 
in line with expectations. 
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Figure 12. Using a critical mass ratio r > 3 to define major mergers does not yield a red sequence, as indicated by colour-magnitude 
diagrams (top row) and luminosity functions (bottom row) for hot gas quenching (left column) and merger quenching with r = 3, 5, and 
10. These simulations were run in boxes of side length 2&h~ 1 Mpc. Even at this smaller volume, hot gas quenching results in an obvious 
red sequence with number densities close to those observed. Merger quenching with values of r > 3 leads to a stronger suppression of 
galaxy growth (fewer bright galaxies, e.g. right column), but does not yield a distinct red sequence with substantial number densities. 
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Figure 13. Relations between quenching energetics and galaxies, for both merger winds (green) and hot gas heating (black). Left: 
Energy associated with each quenching event vs. the stellar mass of the galaxy hosting that event. Since hot gas heating is actually 
continuous, each "event" is amb iguous, but correspon ds to ~ 10 7 yrs. Errorbars show an estimate of the energy output in winds from 
post-starburst galaxies based on lTremonti et ahl j20071 ). Corresponding color-coded histograms of stellar mass are shown at the bottom, 
with blue for optically-selected SDSS AGN. The dotted line shows a scaling E oc M 2 / 3 motivated by the hydrostatic temperature scaling 
with halo mass. Middle: Time-integrated energy from quenching feedback as a function of z = galaxy stellar mass. This is the total 
energy of feedback processes over each galaxy's lifetime, but excludes feedback events that occur in small galaxie s that merge with the 
main progenitor. The dashed line shows the slope of the A/bh — AfsteUar relation for galactic bulges from lBennert et al.l ||2010| ). Right: 
Total ener gy density productio n rate of feedback processes vs. redshift. Squares show a scaled estimate of the total luminosity density of 
AGN from lBarger et al.l i2005l ). 



We track the energy associated with our quench- 
ing implementations on galaxy-by-galaxy and event-by- 
event bases. For hot halo quenching event, we add up 
the change in energy for each gas particle affected, or 
quenched, by the event. With hot gas heating, each heated 
gas particle changes in energy by AE = 3/2OTbAT = 
3/2(Mp ar ticio/(M TO p)) fc sAr, where M par ti c ic is the total 
mass of the particle, fi — .6 is a typical me an molecu- 
lar weight for gas in clusters |Rosati et alj |2002). m p is the 



proton mass, fcs is Boltzmann's constant, and AT is the 
temperature change of the gas particle. Since each heated 
particle is associated with a galaxy from the group finder, 
we attribute its change in energy to the feedback energy 
tally for that galaxy. 

In our merger quenching model, we estimate the feed- 
back energy for each particle as the final kinetic energy: 
E = 0.5M p vl ick . Here i>ki c k is the kick velocity of gas ejected 
in winds associated with mergers. The total feedback en- 
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ergy from a single merger is then summed over all parti- 
cles that experienced the feedback in a galaxy. Thus, for 
our basic model, the total energy in a merger is effectively 



O.SMcoidga^lSOOkms -1 ) 2 , where M, 



coldga 



is the 



total cold gas mass in the merger remnant at the time of 
expulsion. 

In Figure 1131 we compare our quenching energy calcu- 
lations to observations associated with feedback processes. 
In the left panel, we show the energy of each feedback event 
as a function of its host galaxy's stellar mass. Since hot gas 
quenching is actually continuous, there are too many points 
to plot clearly, so we randomly choose 1000 events. The event 
energy rises steeply with stellar mass since more massive ha- 
los have more gas that requires heating to a higher virial 
temperature. Since hydrostatic haloes follow the scaling 
M (x T 3/2 , we overplot a dotted line with E a M^{ Uar (with 
arbitrary normalization). The hot gas quenching events fol- 
low this trend, although with some scatter, suggesting that 
the virial temperature is the most important factor in de- 
termining how much energy is injected in this model. 

Merger quenching shows more scatter and less steep 
dependence on stellar mass, both of which are driven by 
gas fraction trends. The event energy for mergers is solely 
determined by the cold gas mass (since Ukick is the same 
for all galaxies), so the scatter in the relation owes to vari- 
ations in the cold gas mass at any fixed stellar mass. In 
our simulations, more massive galaxies tend to have smaller 
gas fractions, which counteracts the overall increase in mass 
in these systems, leading to a sub-linear relation between 
merger event energy and host galaxy stellar mass. 

For comparison with the merger quenching energetics, 
we estimate the energy output associated with the massive 
outflo ws observed in post-starburst galaxies. iTremonti et al.l 
(|2007l ) measured outflow velocities ~ 1000 km s 1 , and esti- 
mated outflow masses at M w i n< j ~ 10 9_11 Mq. From this we 
estimate the kinetic energy as 0.5M wind v 2 ind , using the full 
range of mass estimates to get the errorba rs, but neglect- 
ing (s ubstantial) variations in wind velocity. ITremonti et al.l 
({20071 ) also measured stellar masses for their sample of post- 
starburst galaxies, and we plot errorbars based on the stan- 
dard deviation of the logarithm of those masses. Our merger 
outflow energies fall at the high end of this observational es- 
timate, entirely consistent given that we used w w ind = 1500 
km s _1 in the simulation. 

We also show histograms of the stellar mass distribu- 
tion of feedback events along the a;-axis. As indicated by the 
black histogram, most hot gas quenching events are associ- 
ated with galaxy stellar masses between 10 10 and IO^Mq. 
Because of the steep cutoff in the galaxy stellar mass func- 
tion, there are few galaxies above this mass, and galaxies 
at lower masses do not have hot halos. In contrast, the 
distribution of stellar masses of merger events (green line) 
is skewed toward the low mass end, reflecting the greater 
frequency of mergers of low-mass galaxies when averaged 
over cosmic time. For comparison we show in blue the dis- 
tribution of stellar masses of AGN host galaxies for AGN 
selected from SPSS using the Baldwin-Phillips- Terlevich 
dBaldwin et al.l Il98ll) diagno stic with SDSS optical spec- 
tra ( Kauffmann et al]|2003al '). The general trend is similar 
though the predicted histogram is shifted to somewhat lower 
masses. This indirectly suggests that AGN activity is indica- 
tive of ongoing hot halo quenching. 



In the middle panel of Figure 1131 we show the time- 
integrated feedback energy as a function of host galaxy z = 
stellar mass. If the feedback energy originated from black 
hole accretion, this enables comparison with final black hole 
masses, modulo corrections for the efficiency at which ac- 
creted mass is converted to output energy, and the coupling 
of that energy to the gas. As with individual event energies, 
the integrated energy rises more steeply with stellar mass for 
hot gas heating than for merger winds. A dashed line shows 
the slope of the Mbh — M B t e iiar relation for local galactic 
bulges from iBennert et al.l |2010T ). assuming that the inte- 
grated feedback energy is tx Mbh - The amplitude is arbitrar- 
ily normalized to coincide with the mean energy for merger 
winds at the low-mass end. The slope for merger winds, 
though slightly shallower than the Mbh — Miliar relation, 
is a good match to the observations. The slope for hot gas 
quenching is also in good agreement at the high mass end, 
but drops away quickly to low masses where little energy 
input has occured. Note that the observed Mbh — M st eiiar 
relation uses the bulge stellar mass, whereas the simulated 
galaxies show the total stellar mass. Hence the dropoff to 
low masses in the hot gas quenching case may partly be 
reconciled by lower bulge-to-disk ratios in sub-L* galaxies. 

Finally, the right plot shows the energy density pro- 
duction rate of quenchi ng as a funct i on of redshift. This is 
analogous to the plot of lMadau et al. ( 1998) of galaxy lumi- 
n osity density evolutio n, and directly comparable to the plot 
of iBarger et al.l (l2005h of AGN luminosity density evolution 
(see also Cowie fc Barged I2004I ). We take data from Figure 
25 of iBarger et al.l for spectroscopicallv identified sources 
with I/2-iokcV > 10 42 ergs s _1 and multiply by their sug- 
gested bolometric correction of 35 to get the total luminosity. 
Then we multiply by a coupling efficiency which determines 
what fraction of an AGN's total luminosity couples to the 
surrounding g as. We choose a somew hat high value of 0.1 
(0.05 is typical |Pi Matteo et al1l2005l ) to match the normal- 
ization of the energy density production rate for the merger 
quenching model, and show the result as squares in the fig- 
ure. Outflows from galaxy mergers in our simulations match 
the shape of measured energetics of luminous AGNs. The 
total energy requirements for hot gas heating appear much 
higher. This may indicate that hot gas quenching in the real 
universe (1) is tied to lower-luminosity AGN, (2) does not 
heat gas up to the virial temperature, and/or (3) is inter- 
mittent with a small duty cycle. 

Overall, our results are broadly consistent with the 
idea that AGN are responsible for providing the energy for 
quenching feedback. Various authors have already shown 
that merger-driven feedback can lead to galaxies that are 
consistent with the M — a relation; our simulations likewise 
appear to do so. However, our simulations also show that 
merger quenching is, in itself, unlikely to be responsible for 
producing red and dead galaxies. The energetic requirements 
of hot gas quenching are quite severe in our current model 
implementation. For instance, such energy injection will ex- 
cessively heat intragroup gas, likely leading to a mismatch 
with observed temperature and entropy profiles; we plan to 
explore this in future work. We speculate that intermittent 
heating, perhaps triggered by an accumulation of cold gas in 
the halo center, might work just as well to quench star for- 
mation with a significantly lower energy budget. Hence while 
forcing gas in hot halos to remain at the virial temperature 
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Figure 14. Quenching rate (black, left axis) and mass growth 
history normalized to z = (green, right axis) of red sequence 
galaxies with M stc n ar > 1O 9,5 M0 in the hot gas quenching model. 
In this model, quenching peaks at z RJ 1, declining at late times. 
About 2/3 of the mass of the red sequence builds up at z < 1, 
while the mass in the blue cloud declines during the same period. 

produces a reasonable red galaxy population, further work 
is required to construct a fully viable physical model. 



5.2 Red sequence build-up over cosmic time 

The build up of the red sequence over time can place strong 
constraints on quenching physics. Here we investigate two 
simple parameterizations: the galaxy number quenching rate 
and the mass growth histories of the red sequence. For both 
quantities, we restrict ourselves to red galaxies with stellar 
masses above 1O 95 M . 

We calculate the quenching rate by comparing each con- 
secutive pair of snapshot outputs from our basic hot gas 
quenching simulation. The quenching rate, dNqu^^ed/dt is 
estimated as AN\ ov , SB { r / At Bna ,p ■ Here, AiVi owss fr is the dif- 
ference in the number of galaxies with specific star for- 
mation rates (SSFRs) below 10~ 12 yr _1 between the two 
snapshots, and At snap is the time between snapshots. The 
red sequence mass growth is the total stellar mass in such 
low-SSFR galaxies, M Ic d(z), normalized to the total stellar 
mass at in these galaxies z = 0, M rc d(z = 0). We also cal- 
culate the mass growth of galaxies with SSFRs above this 
cutoff, Mbiue(z), normalizing to the mass in red galaxies 
M le d(z = 0) to track relative changes in the blue and red 
populations. 

Figure rril shows the results. The quenching rate is small 
at z > 2, after which it rapidly increases. The rate peaks at 
z ftf 1, declining at low redshift. Observations tend to indi- 
cate that the mean formation time of stellar population s in 
massive ellipticals is at z > 2 (e.g. iGraves et al.ll2009l ). In 
our model, while massive galaxies tend to have old stellar 
populations, they do not quench their star formation until 
typically z ~ 1. This may be related to why the color-density 
relation disappears or inverts at z > 1: At low-z, typical cen- 
tral massive galaxies are red, leading to the familiar positive 
color-density correlation, while prior to z ~ 1 central galax- 
ies have typically not quenched and are blue, leading to an 
inverted relation. 

The total mass on the red sequence at z = 1 is about 
1/3 the final value at z = 0; half the mass on the present- 



day red sequence is in place at z = 0.8. This significant 
late buil d-up is slightly m ore than the factor of 2 in stel- 
lar m ass l|Bell et al.ll2004) and galaxy number dFaber et al.l 
120071 ) suggested by observations, although such estimates 
are presented as lower limits. The mass on the blue cloud 
rises at early times, declining at 2 < 1 as many star-forming 
galaxies become quenched. This reflects a decrease in the 
number density of blue galaxies by ~ 20%, substanti ally dif- 
ferent from the empirical analysis of iBlantonl (|2006t ) ■ which 
found no evolution in number density of the observed blue 
cloud to within 10%. By z = 0, the total stellar mass in 
blue galaxies (above our resolution limit of ~ 3 x 10 9 Mq) 
is ~ 70% that in red galaxies. 



6 DISCUSSION 

6.1 Red galaxy morphologies 

In this work we have focused on the transformation of 
galaxy colours, but of course colour and morphology are 
closely correlated in the local universe. Our simulations' 
resolution is poor ly suited to study morphologies (see e.g. 
iBrook et al.l l201ll . for a fuller discussion of why) . In our 
favoured hot gas model, starvation by the hot halo only stops 
star-formation without changing the morphology. Once star- 
formation stops, a galaxy must undergo (possibly minor) 
mergers to disrupt the remaining disk and form an ellip- 
tical. Such dry (gas-poor) mergers may have difficulty re- 
producing the observed phase-space densitie s of stellar or- 
bits i n the centres of some elliptical galaxies (|Hopkins et al.l 
2009), but some fraction of the bulge may have built up 
via dissipational mer gers before the galaxy starved (e.g. 
iBournaud et alj|2007ft . 

One key argument in favor of gas-rich major mergers 
as a transformative process is that they might turn a star- 
forming disk into a quiescent elliptical all at once. How- 
ever, recent observations of quenched late-type galaxies dis- 
favor the merger-only scenario, and suggest that a substan- 
tial fraction of passiv e galaxies must have gone through a 
"pass ive disk" phase dBundv et alj|20icl : Ivan der Wei et al.l 
|2010j). Our results reinforce the theoretical difficulty seen in 
earlier models: after halting star-formation at any given time 
in a galaxy's evolution, something must pre vent future gas 
accretion to keep that galaxy red and dead (|Hopkins et al.l 
|200S| ; ICroton et al.ll2006l ; iBower et al.ll2006l ). 

Finally, we note that our hot gas quenching model does 
not produce passive field galaxies in low-mass haloes (which 
are not dominated by hot gas) , and this may violate observa- 
tional constraints on the relationship between stellar-mas s 
and the fraction of passive galaxies (|Hopkins et al.l r2008). 
Ultimately, morphologies provide an interesting constraint 
on quenching models, but higher-resolution simulations or 
sub-resolution morphological models are required to prop- 
erly apply them. Ins ofar as galaxy struc ture influences star- 
formation histories (M artig et al.l [20091 ). such models may 
be necessary to understand red galaxy formation. 

6.2 The importance of hot halo gas 

Our favoured model relies on the hot gas surrounding mas- 
sive galaxies. This hot corona is a robust prediction of years 
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of theoretical work, yet its existence has not been universally 
established at mass scales below groups. X-ray observations 
have revealed diffuse h ot gas around only some of the most 
massive field galaxies (jMulchaev fc Jeltemal 120101 ) . but re- 
cent simulations predict that the halo gas around ^optical 
galaxies shoul d have X-ray lumin osities too faint for current 
observatories (|Crain et alj|2010al ). The gas may be there, we 
just cannot see it. 

Why should the hot gas be implicated in shutting down 
star formation? The hot gas may be crucial to triggering 
the type of fee dback that best couple s with surrounding gas 
(i.e. radio iets lDi Matteo et al.ll200fj| ). or it may simply be a 
more effective "net" than cold IGM to "catch" the energetic 
output of an AGN. It seems an unlikely coincidence that the 
mass scale where a stable hot halo can form just happens 
to be around the turnover in the stellar mass function, but 
observations do not unambiguously favor a direct connec- 
tion. Among field galaxies with diffuse X-ray gas detections, 
both star-forming a nd passive galaxies are well-represented 
l|Crain et alj|2010bh . This result is not entirely inconsistent 
with our model - after new accretion is quenched, a galaxy 
will remain star-forming for ~ 2 Gyrs while surrounded 
by a substantial hot gaseous halo. Nevertheless, the under- 
abundance of massive star-forming galaxies in our favoured 
model suggests a less direct relationship between hot gas and 
quenching than our simple prescription. Perhaps quenching 
depends jointly on hot gas and the mass of the central black 
hole, which itself is the product of a complicated history 
that may involve galaxy mergers and secular growth. 

6.3 Comparison to other models 

Our on-thc-fly quen ching mode ls bui ld on the post- 
processing models of iGabor et all (2010), and the results 
show some similarities. The color of the red sequence is too 
blue and roughly constant with mass, which we attribute 
to difficulties with supernova metal yields, the redistribu- 
tion of metals by star-formation-driven winds, and/or a 
poor treatment of merger- induced starbursts in our simu- 
lations. We also observe a similar sharp cutoff in the blue 
galaxy mass function in our best model. A major differ- 
ence, however, is that re- fueling of merger remnants prevents 
our merger quenching model f rom yielding a red sequence, 
whereas in IGabor et alj (|20 10T ) we assumed that all future 
star-formation was quenched after a merger. 

With many successes at matching observations, SAMs 
also serve as a useful benchmark for our models. Unlike our 
simulations, typical SAMs already resulted in a substantial 
red sequence even before invoking quench ing feedback mod- 
els (|Croton et al.ll2006l ; [Bower et alll2006h . The red galaxies 
arise due to rapid gas consumption in a starburst associated 
with a merger or disk instability, or due to gas stripping 
as a satellite galaxy falls into a lar ger halo. Our si mula- 
tions do produce quenched satellites l|Dave et al.ll20lil ). but 
probably do not track secular disk evolution properly. Nev- 
ertheless, our merger quenching model shares the generic 
feature with modern SAMs that further ongiong quench- 
ing is required to ensure that the most massive galaxies 
remain red and dead. Different models incorporate differ- 
ent mechanisms for this added quenching, from various ra- 
dio mode feedback pres c riptions to a simpl e halo mass cut- 
off (|Croton et all l200rj ; iBower et al] l200fj ; ICattaneo etafl 



2006; Somerville et all I2008T ) . It is worth recognizing that 
our late-time re-accretion of gas is dominated by recycled 
winds, something that is not dynamically followed in the 
SAMs. 

Cosmological hydrodynamic simulations akin to ours 
have included various feedback processes, but (to our knowl- 
edge) none has previously produced a realistic red se- 
quence and mass/luminosity function. Simulations with 
neither stellar nor black hole feedback produce massive 
galaxies with low star-form ation rates, but t heir stellar 
masses are much too large (K eres et alj [20 09a). Our sim- 
ulations with only stellar feedback can bring the low-mass 
end of the galactic stellar mass function into agreement 
with observat ions, but massive galaxie s form stars at ex- 
cessive rates (|Oppenheimer et al.l 1 2010) . By coupling sub- 
grid Bondi-Hoyle-Lyttleton black hole accretion with mod- 
els for the associated energy output, other simulations 
show promise for matching black hole grow th and sup- 
press i ng star-formation in massive galaxies jSiiacki et alj 

a 



20071; iDi Matteo et al.l [20oi; iBooth fc SchavJ I2009I . |2010| ; 



Degraf et al.ll2010l ; iMcCarthv et al.ll2010l) . but it is unclear 
whether they produce a red sequence as observed. Under- 
standing the differences between these models and ours will 
provide complementary insights into the physical processes 
behind quenching. 

6.4 Star-formation driven winds and quenching 

lOppenheimer et alj (|2010h highlighted the importance of 
wind recycling accretion in our simulations (without quench- 
ing), showing that star-formation in galaxies with M a t e iiar > 
10 10 Mq is fueled mainly by gas that was earlier ejected in 
winds. Because this wind recycling dominates the fuel for 
star-formation in massive galaxies, it is the main fueling 
mode we are trying to quench in this work. 

Our treatment of star-formation driven winds may not 
be physically accurate. Once launched, we decouple the 
winds from hydrodynamic interactions with surrounding gas 
to all ow them to escape the galaxy in a resolution-con verged 
way l|Springel fc Hernou"istll2003l ; ISchave et alj|2010h . Even 
after re-activation of the hydrodynamics, the wind particles 
may not merge seamlessly with surrounding gas. Upon es- 
cape from galaxies, wind particles are cold, dense, and metal- 
rich. While they decrease in density as they move away from 
the galaxy, they remain enriched since we don't include diffu- 
sion of metals among gas particles. Metal-line cooling makes 
them cool faster than surrounding gas, so they stay cold and 
dense, and more likely to return to galaxies. It is unclear how 
realistic the assumed complete lack of mixing is, but we note 
that this wind model has resulted i n agreement with ob- 
servations of IGM metal absorbers (lOppenheimer fc Pavel 
I20061 . 120081 ; lOppenheimer et al.ll2009h . 

Wind recycling, and the possibility of over-recycling, 
has implications for all ou r que nching models. As shown in 
lOppenheimer et all (|2010h and IGabor et al.l i|2010l '). simply 
stopping wind recycling is unlikely to produce a realistic red 
sequence - doing so suppresses star formation too much to 
create enough ~ L* galaxies, while small amounts of cold 
and hot accretion keep too many galaxies forming stars. In 
our merger quenching model, we have emphasized the re- 
fueling of merger remnants after ~ 2 Gyr. Based on the 
origin of star-forming gas explored in lOppenheimer et al.l 
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(2010), most of the accreted fuel comes from wind recycling, 
although a significant sub-dominant component comes from 
cold and hot accretion. Therefore in our hot gas quenching 
model, where we expressly heat all non-star-forming gas, 
most of the gas we must quench had been expelled in winds. 
Hence we caution that the hydrodynamics of winds may 
play a critical role in quantitatively assessing the amount 
and location of quenching feedback needed to produce a red 
sequence. 



6.5 Guidelines for future models 

Given the wide range of physical scales important in galaxy 
evolution, developers of cosmological simulations have little 
choice but to resort to sub-resolution quenching models for 
the foreseeable future. Here we discuss ways of improving 
these models. Any model quenching mechanism requires 3 
steps: 

(i) Criteria to start (and stop) the feedback process. 

(ii) An estimate of the energy budget or luminosity avail- 
able. 

(iii) A method to couple feedback energy to surrounding 
gas. 

The first two steps may be linked in models where the 
feedback energy varies over time: e.g. the black hole ac- 
cretion rate varies by several orders of magnitude between 
states where the feedback is important and negligible. In our 
current hot gas quenching model, (i) is the existence of a hot 
gas-dominated halo, (ii) is effectively unlimited, and for (iii) 
we heat all the non-star-forming gas in the FOF group. We 
have seen that this successfully creates a red sequence of 
galaxies, but that it leads to a dearth of star-forming galax- 
ies with Msteiiar > 10 10 ' 5 Mq and requires too much energy 
input. Below we discuss each of these steps in the context of 
radio mode AGN feedback. 



6.5.1 (i) Triggering mechanisms 

Triggering radio mode feedback requires accretion onto a 
supermassive black hole. In massive (cluster) haloes, the 
black hole ma y accrete directly from th e hot halo in a 
Bondi process (Qu ataert fe Naravan|[2000h . or a black hole 
may be f ueled by col d clouds that condense f rom the hot 
medium (|Sokeij I2OO6I : iMcNamara et al.1 |2010|) . Bondi ac- 
cretion of hot gas may be more conducive to the forma- 
tion of a power ful radio jet than rap id, radiatively effi- 
cient accretion (|Di Matteo et all 12000). The formation of 
a jet may additionally req uire significant black hole spin 
l|Naravan fc Qu ataert 2005;) . Based on the present work, the 
presence of hot gas seems to be crucial for triggering radio 
mode feedback. 

Models have used a variety of triggers for AGN feed- 
back. SAMs typically assume that feedback operates in 
hot, massive haloes. Hydrodynamic models usually adopt 
a Bondi-Hoyle prescription, or a simple scaling with gas 
mass in the black hole vicinit y (|Vernaleo fc Revnolds!l2006l ; 
Brug gen fc Scannapiecoll2009l ). to control black hole accre- 
tion and thus triggering. 

Observationally, most X-ray bright clusters, groups, and 
elliptical galaxies host some level of central radio AGN, with 



AGN more likely in systems with shorter central cooling 
times, and more luminous in more massive systems l|Burns! 
ll99d : lMittal et al.ll2009l : ISunll2009l ; |punn et al.ll201ol ). Out- 
bursts that induce X-ray cavities typ ically occur every 10 7 ~ 8 
years with e nergi e s ~ 10 55 - 60 ergs JVoit fc Donahue! 120051 : 
iDavid et al. I l2010l : iRandall et alj|2010l ). These observations 
suggest a triggering mechanism linked to gas cooling in the 
central regions of a hot halo, which may be a resolvable 
triggering mechanism in cosmological simulations. This may 
provide a more observationally-calibratable method for trig- 
gering AGN feedback than Bondi-Hoyle accretion. 



6.5.2 (ii) Available energy for feedback 

The energy of feedback mainly depends on the black hole 
accretion rate. Most models assume some fraction of the 
accretion rest-mass-energy acts on surrounding gas. The ac- 
cretion rate is typically determined through simple scalings 
with black hole and halo mass in SAMs and through Bondi- 
Hoyle models in simulations. 

Observations suggest kinetic power output from central 



bursts of 



10 45 erg s 1 over 



10 years in some clus- 



ters, consistent with total energy in inflated bubbles (~ 10 
ergs). As mentioned above, radio emission scales with clus- 
ter size, but it is not clear to what extent the burst energies 
depend on halo properties. The fact that few very massive 
galaxies have enough star-formation to give them blue colors 
suggests that there is enough heating pow er to offset cooling. 
This may imply a self-regulating cy cle (cf. lCiotti fc Ostrikerl 
1 19971 . l200ll : IChurazov et all [2005). Hydrodynamic models 
that directly track black hole accretion or use some esti- 
mate based on the mass within the central regions of the halo 
can a chieve this self-regulation fe.g. lBriiggen fc Scannapiecol 
120091) . 

In our hot gas quenching model, the required energy 
budget is large, although it is not clear that significantly 
less energy input could be equally effective. It is also difficult 
to exceed the theoretical maximum energy output of black 
holes; models typically assume efficiency of mass-energy cou- 
pling to surrounding gas of less than 1%. So the energy bud- 
get remains an interesting constraint, but the most attrac- 
tive scenario is that black holes emit just enough energy to 
remove or starve their cold gas reservoir, at which point they 
become quiescent again. 



6.5.3 (iii) Coupling of feedback to gas 

How AGN energy couples with the surrounding medium 
is another matter of debate. Supe rsonic jets induce bow 
shocks in th e surr ounding gas (IVoit fc Donahue! 120051 : 
IRandall et al] |2010| ). and subsonic jets entrain mate- 
rial and therma lize their energy via turbulent mixing 
(|De Yound |2010D . The outflows can apparently evac- 
uate ~ 10 kpc-sc ale "bubbles" in hot intragroup gas 
(Birz arTet al.l [2004), and the bubbles then share their 
energy at scales ~ 10 — 100 kpc via buoy a nt uplift, mixing 



and sou nd waves dChurazov et al 
iNulsen et all 120071: 



2001; Voit fc Donahue 



Scannapicco & Bruggen 



1 20081: iBriiggen fc Scannapiecol 120091 : iFabian et all 120031 : 
Zanders fc Fabianll2007l , l2008h . 



While idealized hydrodynamic models have directly 
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examined the i nteraction of jet s and bubbles with 
the I CM (e.g. iRuszkowski et alj 12004 iBriiggen et all 



20051: iHeinz et all 120061: IScannapieco fc Briiggenl l200Sl: 
Briiggen et all 120*091 ; iMorsonv et al.ll2010l ; iMendygral et all 
20 Ilk SAMstypicallv assume that feedback directly offsets 



the cooling for simplicity. Cosmological hydrodynamic sim- 
ulations usually inject the feedback energy into one or a few 
resolution elements near the black hole, relying on the hy- 
drodynamic calculatio ns to distribute the energy to large 
scale s in the halo (e.g. iBooth fc Schavell2009l ; iTevssier et all 
2011). In our models, suppressing star-formation in accor- 
dance with observations requires heating gas throughout the 
circum-galactic region, on scales of > 10 kpc, to prevent the 
formation of cold clumps. 

In summary, a path for modeling quenching feedback 
would be as follows: 1) Trigger feedback when the central 
cooling time drops below e.g. 1 Gyr, or some fraction of 
the hubble time; 2) Calculate energy for feedback as some 
fraction of the rest-mass energy of cold or cooling gas in 
the central regions; 3) Increase the entropy of the gas with 
some radial dependence in a manner consistent with higher- 
resolution simulations. We will explore models based on 
these principles in future work. 



7 SUMMARY AND CONCLUSION 

With the goal of building a realistic red sequence of galax- 
ies in cosmological hydrodynamic simulations, we have im- 
plemented novel mechanisms for quenching star-formation 
in the simulation code GADGET-2. By identifying merg- 
ers and halos on-the-fly within simulations, we implement 
and test various models for quenching feedback related to 
these processes. While our work is motivated by feedback 
resulting from AGN, we explicitly avoid examining black 
hole growth in order to concentrate on constraints from the 
massive galaxy population. 

Motivated by observations of massive outflows from 
post-starburst galaxies attributed to quasar feedback, we 
implement 1500 km s _1 superwinds in the remnants of 
galaxy mergers that expels all the cold gas. To do so, we 
first use an fof group finder to identify galaxy mergers on- 
the-fly, and then give a velocity kick to all the star-forming 
gas. We have shown that this quenching mechanism alone 
does not produce a red sequence in our simulations. Even 
after all the gas is expelled from a merger remnant, new sup- 
plies of gas accrete from the IGM to re-fuel star-formation, 
typically within ~ 2 Gyr. 

Motivated by observations of radio AGN and X-ray cav- 
ities in the hot gas of galaxy groups, we add thermal heat- 
ing to hot gas in massive dark matter haloes. We calcu- 
late the hot gas fraction /hot in the halo, and if /hot > 0.6 
(roughly corresponding to halo masses > 10 12 Mq), then we 
heat all its gas outside of the ISM to the halo virial tempera- 
ture. By keeping the surrounding gas hot, we starve galaxies 
embedded in hot haloes of new fuel for star-formation. We 
have shown that this hot gas quenching model yields a red 
sequence whose luminosity function provides an excellent 
match with observations of local galaxies (Figure [5]). 

Our main results are: 

• Galactic-scale outflows triggered by mergers (i.e. quasar 
mode feedback) do not produce a substantial red sequence 



on their own because gas accretion from the IGM re-fuels 
star-formation within ~ 1 — 2 Gyr. 

• Adding thermal energy to hot X-ray gas around mas- 
sive galaxies (analogous to radio mode AGN feedback) suc- 
cessfully produces a red sequence whose luminosity function 
matches observations. 

• This heating must occur around satellite galaxies em- 
bedded in the hot gas of their parent haloes to match the 
faint-end red galaxy luminosity function. 

• A combination of hot gas heating with merger-triggered 
outflows may be empirically motivated and perform as well 
as the heating-only model, but hot gas heating is the crucial 
required component. 

• Our simple hot gas heating model produces somewhat 
too few massive blue galaxies, possibly owing to the sharp 
truncation in gas accretion onto galaxies as soon as their 
haloes are dominated by hot gas. 

• As in lGabor et al. I fcoich , our baseline model produces 
a red sequence that is too blue and too shallow, likely ow- 
ing to issues related to enrichment. We empirically recali- 
brate the metallicities to obtain agreement, which requires 
a metallicity increase up to x2 in the most massive systems. 

The overall success of reproducing the observed red se- 
quence and associated luminosity function is a first for cos- 
mological hydrodynamic simulations. Doing so has already 
elucidated stringent constraints on how quenching must op- 
erate under the scenarios explored. While hot halo quench- 
ing appears to be necessary and sufficient to reproduce ob- 
servations of red galaxies as well as any model at z = 0, our 
current simplistic implementation requires more energy than 
is thought to be available and likely overheats surrounding 
gas. We have yet to explore details of redshift evolution, 
clustering, or the impact of quenching on the surrounding 
intergalactic gas; these may motivate variants on this simple 
scenario. The beauty of implementing quenching models self- 
consistently within large-scale hydrodynamics simulations of 
galaxy formation is that such models open up a host of new 
avenues to constrain quenching physics. 

In future work, we aim to develop more physically con- 
sistent models for quenching star-formation. To do so, we 
advocate connecting well-resolved structures and processes 
in the simulations to feedback processes that operate be- 
low the resolution scale. Such sub-resolution models will re- 
quire guidance from observations and higher-resolution sim- 
ulations of individual galaxies or clusters. We outline a par- 
ticular path forward for this based on current intuition that 
we will explore in future work. By combining insights from 
advancing theoretical work on all scales and advancing ob- 
servations across cosmic time, we hope to continue refining 
our models to better understand how massive red and dead 
galaxies come to be. 
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